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¢ : Edit History for Version 1 of MTHSHEXP 

9000 54; 1-90) - Adapted from MTHSGEXP version 1-002. JAW ea 

0000 55 ; 1-002 = Added Large argument logic to prevent lose of s or) ficence 

0000 56 ; for arguments with absolute value greater than 2**7 

0000 af 3 RNH 24=JUN=81 
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2090 23 -SBTTL DECLARATIONS ; Declarative Part of Module 
00 61; 
4 $¢ ; INCLUDE FILES: MTHJACKET.MAR 
it be 3; EXTERNAL SYMBOLS: 
itd 66 -DSABL GBL ; Declare all EXTRNs explicitly 
0000 67 -EXTRN MTHSSSIGNAL ; SIGNAL SEVERE error 
44 rt: -EXTRN MTHSSJACKET_TST 3 Test to see if called with CALL or 
9900 70 .EXTRN MTHSK_FLOUNDMAT : Underflow error code 
B00 4 -EXTRN MTHSK_FLOOVEMAT ; Overflow error code 
4h 73 ; EQUATED SYMBOLS: 
0000407C 0000 75 ACMASK = “M<IV, R2,_R3, R4, R5, RO> ; «ENTRY mask + int ovf enable 
00007D12 44 £6 X_16LOG2E = *x7D12 ; Extension for operand in EMODH 
0000 78 ; MACROS: 
0000 79 ; 
it + SSFDEF : Define SF$ (stack frame) symbols 
4 6 3 PSECT DECLARATIONS: 
00000000 84 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
0000 85 ; Program section for math routines 
0000 86 ; 
0000 87 ; OWN STORAGE: none 
0000 88 ; 
0000 89 ; CONSTANTS: 
0000 90 
0000 91; 
0000 38 ; Table of 2**(1/16) for I = 0 to 15 
0000 «94 
00000000 00004001 0000 95 TABHI: .LONG %*X00004001,*x00000000 =; 2#*(0/16) = 
00000000 00000000 0008 96 eLONG °%X00000000,°x00000000 =; _0.1000000000000000000000000000000000E +01 
989086CF 08554001 oats 97 LONG Ree er Ht LP 3: 2e2(1/16) = 
42AAB718 sozel Ose 001 98 LONG Rs 334: A os 14h : 0.104427378242741 38403219664 78739929E +01 
D51783C7 17284001 0020 Vo) LONG *X172B4001,°XD51783C7 ; 2*(2/16) = 
BSAC 308 F7C8ADCD Bap 100 LONG Bore m8 30C206 : _0.1090507732665257659207010655760707E +01 
5623A6E7 23874001 0030 £101 LONG *%X23874001,°X5623A6E 3 202(3/16) = 
5CB6B1C1 AD eEee iF 19 LONG “X1FAD8 § 2x3 CRG) CI : _0.1138788634756691653703830283841511E+01 
187 te ORs 4 of 01 040 10 LONG mF 1,°X1B71E0A3 3 2%2(4/16) = 
5€864630 8D5A52DE Baek 104 LONG “X8D5A52DE,°xX5 reread 3: _0.1189207115002721066717499970560475E +01 
Se rhe eet bey 105 LONG at 7 EAS Ae 4 gee 3 2%%(5/16) = 
D774 at 41222358 005 106 LONG %*X4122235B8,°XD7745E1 3: 0.1241857812073484048593677468726595E +01 
62A2AD5 eed: bce 107 LONG mes ad 2A2AD5 3 2%*(6/16) = 
sb ireCe ; 7A71D4 38 108 LONG *X397A71D4,°xX Be IDE Ce ; 0. Ger pF AGS 1OOMES TESST IAT TTT IRIE TON 
5427DD4 aitttte 0 109 LONG “*X nt ett tare 427004 ss Qee( £183 = 
Pcresy 9301958C 0078 110 LONG *X930195 Te rE be 3 eat Fr §5546936892728298014740140702E +01 
f gf 67 6A094001 it 111 LONG “*X6A094001,°XF act 67 ; 2%(8/16) = 
£A951 ger oC? 8 0088 ie LONG “*X er OC? 8, “XEA 96 : 0.1414213562373095048801688724209697E +01 
B018475E 7A114001 64 11 LONG °*X7A114001,°XB018473E ; 2%#(9/16) = 
ita; 5102607D 009 114 LONG “*xX51026D Re At : 0.147682614593949931 1 386907480374049E +01 
AA0D5422 BACE4001 OOAO 115 LONG “*X8ACE4001,*xXAA0D54 : 2e#(10/16) = 
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C9BBA192 7CSS5SB5BA AB 11 -LONG “*X7CS5S5BSBA,*XC9BBA192 ; 0.1542210825407940823612291862090734E+01 
aa 9€494001 Oe0 113 - LONG “196290001 “ASE O9T BZA 5 awe erat? = 
BE60 F €46B01C7 B 118 LONG Rss 14 Oe BE6071F ; _0.1610490331949254308179520667357400E +01 
AD A338 ees} 3" 119 LONG “XA 0 1,°XAD A998 3 2ee(4 119) z 
gysalr? 754DD5E8 00C 120 LONG “*X734DD oo akee Al : _0.1681792830507429086062250952466429E +01 
29 43 bi zpey 1 $08 121 LONG *XC1994001,°x5529 43: 3 2ee(4 At = 
1BA62A0 galt 2 00D 1 § LONG “XOCBIC BN 9 AO s .0.17 $¢ 2160373299483112160619375313E+01 
BA4BBDCF 05814 a QOEO 1 -LONG “*XD5814 Ue a DCF ; 2%#(14/16) = 
EODDEB66 AOSA725D OOE8 124 ONG “XAO05A7 Boe OBOE RSS : _0.1834008086409342463487083189588288E +01 
49S0DFA2ZA EA4A4001 bore 125 LONG “XEA4A4 bas ta A 3 2%%(15/16) = 
DB3018F5 F73A9858 0100 : $ LONG “XF73A9858,*°XDB3018F5 ; 0.1915206561397147293872611270295830E +01 
100 128; 
0100 129 ; Each entry in TABLO gives the difference between the corresponding 
0100 130; entry in TABHI and the power of 2 which it approximates, i.e., 
3198 132 ; TABHICn] + TABLOCn] = 2**(n/16), to 226 bit accuracy. 
100 133° 
00000000 00000000 8180 134 TABLO: .LONG %*x00000000,*x00000000 =; 0.000000000000000000000000000000000E +00 
00000000 00000000 0108 135 -LONG *x00000000,*x00000000 
DD3O0SBAF got ere 0110 136 . LONG eR EE + -9.374520292280427421957567419731049E=-35 
64EB15EA 7 (36530 0118 137 LONG “X7 te Do AGSEB I 2EA 
7FDFD542 E4803F8F 0120 138 LONG “XEG8O3FBF,“°X7FDFD542 =; 9.112493410125022978511686101672356E-35 
AF748F09 A4FC5D61 oI $5 139 LONG “XA4FC5D61,“XAF748F 09 
65BF35EA BI3SFBFSF 0130 140 LONG “XBI3FBFSF,“X65BF35EA =; -8.148468844525851137325691 767487803E-35 
1A8576FE pbetg 0138 141 LONG Mb ets tps ahaa 
BC9N3DB8C 2134BF8D 0140 145 LONG “%X2134BF Bo aCe 20 et 3 -1.359830974688816973749875638246305E-35 
46820E93 O05 GOA 0148 14 LONG “XDOES60AF ,*xX46820E9 
F6OEFIF51 174DBF8B 0150 144 LONG “X174DBF8B,*XF6EFIF51 ; =-3.283170523176998601615065965334027E-36 
00C20344 AEEF660E B28 145 LONG ameecreede*aneerey aad 
EA6345D01 FCOCBF8D 0160 146 LONG “XFCOCBF8D,*XEA6345D1 =; -2.391474797689109171622834301602562E-35 
CEBO4EC6 7D1BA860 0168 £147 LONG A eats TiS t, 
4A035F20 76233F8F 0170 148 LONG “X76233F8F,°X4A055F20 ; 7.036756889073265042421737190671412E-35 
AOC8B2F EF6D81A7 0178 # 149 LONG one rehtres ‘antentente 
O5D4EB7B F4F83F8F 0180 150 LONG “XFGFE3F8F,*XO5D4EB7B =; 9.422242548621832065692116736394105E-35 
2C68BD 42674900 0188 #151 LONG *%X426749DD,*°X642C68B8D 
SDESAD9A 7BDO3F8F 0190 13 LONG “X7BDO3F8F,*XSDESADIA ; 7.143528991563300614523273615092530E-35 
5F016CD2 mm 14C 0198 15 LONG Re EL shai ecDg 
| 3 OFS F3763F ‘ Q1A0 3 =: 154 LONG “XF376 RF. 3999B0F : 4.696933478358115495309739213201925E-35 
Doeg fone A36C825 af 155 LONG “%*XA36C8251,°x0 f 686 
F7A386BC C894BF8F 0180 156 LONG “XC894BF8F,°XF7A586BC ; -8.587318774298247068868655935103793E-35 
33A8A216 BD2E9E4D p188 157 LONG “XBD2 Ze ade XSaARAC IS 
D7A9EB99 FEFOSF8F O01C0 158 LONG ontee 3FBF,°XD7A9EB99 3869; «=69..609733932128012784507558697141504E-35 
BEDoA 92 58B5102C git 159 LONG “x5 ? 10 Sons DOA 93 
091259 SESI3SFBF 01D 160 ONG “*X9E513F8F,°X28D91259 ; 7.792430785695864249456461125169745E-35 
Saree CC 00109849 01D 161 LONG “XxX ag 3 F9ESCC 
8FDD82A $708 OF 3 196 LONG ane ¢ BFSF ,*x anes A: -6.472995147913347230055214575612342E-35 
70BA4139 4TH 1€ 16 LONG “XxX fan ’ons 4139 
610CEA20 F86D3F 8D 4 164 LONG “XF86D5F8D,°X61 CEA 0 3; 2.371815422825174835691651228302996E-35 
65278800 3EB/7EEE1 irs 192 LONG “X3EB7EEE1,*°X65278800 
at 
09 4 : Constants used in evaluation of polynomials - small arguments 
0 171° 
DECiF87C AESS3FD8 8 172 HXPTB1: .LONG “KXAESS3FD8,“XDECTF87C 8; 7.6476304910593727523294512938331 16E-15 
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48944388 DCSD9BEF 0 173 -LONG “XDCSD9BEF,“X48944388 
75at0496 BS3P ay 1 174 . LONG “1953 a} sH7BATOG 6 3; 1.147149241385367053304343571061069E=11 
2 §pBA3D SAF 195 8 1 175 eLONG “X9AF19596,*X2 ¢bBA D 
€138612E ghee a3 176 -LONG “%X6124 a3 o"XC1358612E =; 1.605904380491393922264232745268590E-10 
1GOEFFAF 304857 F 17 ONG %*x304857 Featge jh 
SFSIBEFS yEED re 178 LONG “*X EFD FE4,“XOFSIBEFS ; 2.087675694046229638056530752087440E-9 
gp A1B98 3FOF3F 3 17 LONG “X3FOF Dt A1B98 
ageqert aces FE 4 180 » LONG Rly FE7,*Xx ag3e6cF : 2.505210838544287614604977125833492E-8 
it co? 5A7615B2 4 181 ONG “X5A 01288 xRe 5 co? 
AO8FB77 oreo secs 5 1 ¢ LONG “X27E4 Cee aneeee ng 7 =; 2.755731922398745755905032035036362E-7 
perro DC546DE 5 18 LONG aN TOO RODE. mpecee 8 
6C733A 4h ad $9 Is » LONG anges BER aioe 33A5 : 2.755731922398589062864821927508175E-6 
EDA14A28 8 8E 6 185 ONG “X8C8158E3,*XEDA14A28 
1A0101A0 AQIASFF1 4 186 LONG “XAQIA3FF1,*X1A0101A0 ; 2.480158730158730155861058495567987E-5 
C141B47C S2309FF7 027 187 LONG “X52309FF7,°XC141B47C 
1A0101A0 AOIA3SFF4 0280 188 -LONG “XA saree a antesannad 3 1.9841269841269841269841550026571 36E=-4 
339AD98C O1BBAOIA 0288 # 189 «LONG “X01BBA01A,“X339AD98C 
190) C166 6C163FF7 8 90 190 LONG mFS At Ow BERL 8 3; 1.38888888888888888888891 74 16932215E-3 
431467A7 3h dB Ht 191 LONG Fh dD 1467A 
11111117 11113FFA O2A0 135 LONG “X117713FFA,°X11111111 «= s_- 8. 333333333333333333333333331592122E-3 
S800 0EE 11111111 O2A8 19 LONG ada ROB AO TOES 
55555555 a2eeeec’ 0280 194 LONG aneeeeaee ge ak 555555 = s_-_ «4. 166666666666666666666666665273693E-2 
Pett Thy 55555555 0288 195 LONG *x5555555 ange h 2238 
5555555 ry aa | 0200 196 ~ LONG angeeeatc esate? 5555 3 1.666666666666666666666666666666671E-1 
556A5555 55555555 02C8 197 LONG %*X55555555,°x556A5555 
00000000 00004000 0200 198 LONG *%*X00004000,*xX00000000 =; 5.000000000000000000000000000000025E-1 
00180000 00000000 0208 199 LONG “%*Xx00000000,*x001B0000 
00000000 00004001 02E0 200 LONG %*xX00004001,*x00000000 =; 9.999999999999999999999999999999999E-1 
00000000 00000000 O2E8 201 LONG *%*x00000000,*x00000000 
000000 00004001 O02F0 soe LONG %*X00004001,*x00000000 =; 9.999999999999999999999999999999999E-1 
00000000 00000000 O2F8 20 -LONG *%*xX00000000,*x00000000 : 
0000 0300 Soe HXPLN1=<.-HXPTB1>/16 3 No. of entries in table 
0300 05 
0300 206; 
0300 207 ; Constants used in evaluation of polynomial - regular args 
0300 208 ; 
Bs 
16B5D6A9 314B3F9D itd 11 HXPTAB: .LONG “X314B3F9D.*X16B5D6A9 ; 9.407666858107355479744879542837367E-31 
F4E 26068 gece FD3 b3i6 \¢ -LONG “X8EE63FD3,“°XF4E26068 
0891E189 81643FA5 3 $} eLONG “*X81643FA5,°X0891E189 ; 3.040223120834586251881679850732127E-28 
BB46C345 1 bas Hh oa! 14 LONG “X1B98895D,°XBB46C345 
ornater ¢ tts 308 1? thee oa eer ce atk AOSSRr : 9.12281307564508424034 76385 36340806E-26 
z604e 780 GBtasras 35 2 ae “HEBCASTES = 6c? g > 2.5269959389687320708599625 78558151E-23 
Bberst 38 g4cr Seep 340 4? LONG - Ar tia : 6.41640474059492391242500600785501 7E-21 
5 Beans as 255r¢5 350 1 LONG “XB5 53 c$: x3 80 D3 : 1.481106447934839518765528903490102E-18 
D71F0E10 Fae 789 5 ¢ LONG *XEA4B07 eeane 1F0E10 
4 $3 A ef 44) 6 LONG “XxX ef 4 Oe Se A: 3.076970295548511084073493860214445E-16 
41D8C305 6C8735FD ° 4 LONG “*X6C875FD7,°xX41D8C505 Ni 
pets, rete 4} 4 5 Lowe md Be? “HOF ¢AGEBO : 5.682086126528620823749671 267689669E-14 
ec? i «309 Ay as ; LONG “Ha809 SFDC 16¢ a : 9.181219573844438764117422957672360E-12 
A6 eS ROBT SFES 0 9 LONG “Xx Dar Sres “XAOS Pep : 1.2715871950558131621753914 36085248E-9 
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IEESET IC 17F71107 9 0 «LONG “X17F7 1OT ox SESE TOC 
ASE 7B6F ZBOASE EA A 1 . LONG a th FEA, “XA4E/7B6F 3 1.467610032291942926327372981814470E-7 
87FF 83 coer 2 A ¢ LONG CBB fe »“X87FFCBS 
4A0B8D § 3FF 8 LONG “HC6803 FF Fees D7 3 1.355080777949745604324762299029330E=5 
4F12BB2 3A76F 8B B 4 LONG “*X3A76F8B3,°X34F 12BB 
CaBOrE EBFB3FF6 C 5 LONG one OC ESt DO Nes SDFF 3 9.383847928089871575529346217317675E=4 
73CV5EC a sep Gaze C 8 LONG “XOQEGSEF o "X973606EC 
ASO SEE Sees 443 D LONG Xeeee Ee AKAs 4 4+ : 4.332169878499658183857700759113603E-2 
07E66730 95C7F35 D 8 -LONG “X93C7F35 m.t'443 730 
00000000 00000000 E 39 LONG *x00000000 , *x 0000000 ; 0.000000000000000000000000000000000E0 
00000000 00000000 O3€ 40 -LONG *xX00000000,*x00000000 
0000000F 8 : ‘] HXPLN=< CHXPTABD/16 :; No. of entries in table 
O3F0 | 528 
O3FO 44 H_16L0G2 Fs ; LOG2(E) * 16 
B82F7652 71544005 tia 45 ONG *X71544005,*XB82F 7652 
D23AFDAO 7D0FE177 48 $78 "TONG “X7D0FE177.*XD2 3AF DAO 
0400 248 H_LN2_0V_16_HI: : 
00006730 93C7F357 A3S9E2FEF 62E43FFC 84 sep SOCTA *X0000673093C7F357A39E2FEFOZE43FFC ; Hi 95 bits of Ln2/16 
0410 $2) H_LN2_O0V_16_L0: 
069E16C5 4C5B9339 79A157A0 F97B3F96 brie $2 SOCTA *X069E16C54C5B933979A157A0F97B3F96 ; Low bits of Ln2/16 


H 16 
MTHSHEXP ;_H Floating Exponential Function 16-SEP-1984 01:35: AX/VMS Macro V04-00 P 7 
1-002 MTHSHEXP = Standard H Floating EXP ~SEP=-1984 94:32:88 EMTMRTL. SREIMTHHEXP. MAR: 1 ~~ (4) 
? 8 2? -SBTTL MTHSHEXP = Standard H Floating EXP 
hs 
158 8 ;++ 
“ 23 ; FUNCTIONAL DESCRIPTION: 
r 6 1 2 EXP =H floating point function 
r e§ ; Uses a Chebyshev approximation, with overhang on Last step. 
04 65 ; 
8 8 $6 ; CALLING SEQUENCE: 
0420 68 : exponential.wh.v = MTHSHEXP(x.rh.r) 
0420 69 ; 
0420 70 ; -or- 
0420 71; 
8 5 2 ; CALL MTHSHEXP(exponential.wh.r, x.rh.r) 
0420 274 ; Because an H-floating result cannot be expressed in 64 bits, it is 
0420 3 3 returned as the first argument, with the input parameter displaced 
0420 76 ; to the second argument, im accordance with the Procedure Calling 
0420 tt 3 Standard. 
0420 78 ; 
Be 8 ; INPUT PARAMETERS: 
00000004 0420 81 LONG = 4 :; Define sengierd multiplier 
00000008 be50 Hf x = 2 * LONG ; Contents of x is the argument 
0420 84 ; IMPLICIT INPUTS: none 
0420 85 ; 
be 8 $o6 3; OUTPUT PARAMETERS: 
00000004 Se 3 288 exp = 1 * LONG 3; Contents of exp is the result 
93 44 3 VALUE: H floating exponential of the argument 
0420 336 > IMPLICIT OUTPUTS: none 
0420 93 ; 
Bg 3¢ 3; SIDE EFFECTS: 
0420 296 ; Signals: MTHS_FLOOVEMAT if X > 11355.83 with reserved operand in RO/R3 
0420 97 ; (copied to the signal mechanism vector CHFS$L_MCH_RO/R1 Y LIBSSIGNAL). 
4 8 98 ; Associated message vy “FLOATING OVERFLOW IN-MATA LIBRARY’. Result is 
4 99 ; reserved coereng -0.0 unless a user supplied (or any) error handler 
420 0 ; changes CHFSL_MCH RO/R1. 
420 13 MTHS_FLOUNDMAT if X =< -11356.52 and caller has hardware enable set. 
420 ¢ 3; The result is set to +0.0. Associated message is: ‘FLOATING UNDERFLOW 
Be 8 i i IN MATH LIBRARY" 
0420 5 ; NOTE: This procedure disables floating point underflow, enable integer 
04 88 ; overflow, causes no floating overflow or other arithmetic traps, an 
" o3 3; preserves enables across the call. 
4 9 te- 
4 10 
4 11 
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407C¢ 18 sENTRY MTHSHEXP, ACMASK 3; Standard call-by-reference entry 
1 ; Disable DV (and FU), enable Iv 
14 MTHSFLAG_JACKET ; Flag that this is a jacket procedure 


6D, 00000000°GF  9€ MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 


handler 


in case of an error in special JSB 
routine 

RO/R3 = user’ arg 

RO/R3 = epectet HEXP(RO/R3) 

Store result in firat argumen t 
Return to caller 


315 
415 

17 MOVH ax(A RO 
318 BSBB MTHSHEXP. R6 
33 MOVO RO, @exp<AP) 
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1-002 Reo= Special HEXP rout ine oreep= 138 31:30:85 EMTANTS eacdaruincxoomar:1 29% (2, 
~-SBTTL MTHSHEXP_R6 = Special HEXP routine 
: Special HEXP = used by the standard, and direct interfaces. 


: CALLING SEQUENCE: 
: save anything needed in RO:R6 


=z 
“4 
=z==z 
bad 
zn" 
m— 
xo 
ve 


MOVH ans 3 Input in RO/R3 
JSB MTHSHEXP_R6 
return with result in RO/R3 
; Note: This routine is written to avoid causing ony integer overflows, 


: floating overflows, or floating underflows or divide by 0 conditions, 
; whether enabled or not. 


; REGISTERS USED: 
RO/RS = H=-floating argument, then H-floating result 


RTHSHEXP_R6: : ; Special HEXP routine 


The preliminary test for overflow works as follows: First, the sign is 
taken away, leaving just the biased exponent. Then, 16384-4 (bias-4) is 
subtracted, peav ing an exponent biased by 4. This rebiased exponent is 
compared against 18. The comparison can have 3 outcomes. If the rebiased 
exponent is now negative, this means that the true expcnent is < -4 = this 
is a BLSSU test. If the rebiased exponent is pos) tiye- but greater than 
18, then the actual value is greater than 16384, which is guaranteed 
overflow or underflow, depending on the sign of X - this is a BLSS test. 
Otherwise, X is somewhere in the range for the standard evaluation, and 
flow continues. 


else, -4 < unbiased exp < 15 

no exceptions in EMODH or APPROX 

CMPW R4, #*xX400C Check for loss of significance in 
EMOD ( {Xi >= 2##11) 


BLSS EVAL ; No loss of significance 


54 50 8000 8F AB Bicw3 #*x8000, RO, R4 ; Preliminary test for over/under flow 
: R& = exponent bits only 
55 54 3FFC 8F A3 SUBW3 #*X3FFC, R4, RS ; RS = 4 + unbiased guponent 
he i. CMPW #18, R5 ; Unsigned compare with 18 
71 1F BLSSU SMTST ; To more tests if LSSU 


400C 8F 54 861 


iXi >= 2ee4, EMOD will lose significance so the interger and fractional 
parts of X*16/ln2 must be obtained in seperate steps. 


SSS223F5. Eek lk a ak tal ae a ee aE ahah ah teh teh ah ach ah ah ahah ahah ahah ahah ahah ahve ab as tah ahah ah ah ab ah as abias asa! 
FAD AVE WLM MM MMM OO NINN NDI ODA AAA AAA AAA AAA AAA AAA AAA AA AAA A AAA 
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F 65FD MULH3 H_16L0G2_E, RO, -(SP) ; Get_integer part of X*16/\n2 in 
56 gf 6AFD CVTHL (SP), R6 ; R6 T3) as a longword and in 
6E 56 6EFD CVTLH RO, (SP) : (SP) in H format 
; Get fraction part of X*16/\n2 = 
7E 6€ AO AF 65FD MULH3 H_LN2_OV_16_HI, (SP), -(5P) 
50 8E 62FD SUBH2 (SP) +> RO >: 16/Ln2el X = (14d) *Ln2/16 J 
6E AG AF 64FD MULHO H_LN2_OV_16_L0, (SP) ; — in RO/R1. 
50. BE 62FD SUBHe }§=- (SP) +> RO ; 
50 FF7C CF 64FD MULH2 H_16L0G2_E, RO ; 
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1-002 MTH EXP. RO"= Special HEXP routine 6-S 71986 9 a ¢§ YRTHRTL gh TATHHEXP.. MAR; 1 . 
0B sil are 340 BRB APPROX 
56 50 7012 8F FF74 gf 74FD anf 5 EVAL: EMODH 4H_16L0G2_E, #X_16L0G2E, am R6, RO 
4 ¢ : Get X*16*LOG2(E) with 
4 : integer part in RO (=161+J) 
be : : fraction in RO/R3 
FE77 CF OE& 50 75FD Bf $ APPROX: POLYH RO,#HXPLN=-1,HXPTAB ; Use Chebyshev series 
4 ; with last coefficient 0 
8 4 8 ; so that last ADDH has overhang 
54 56 FFFFFFFO 8F CB 82 9 3 BICL3 #16, R6, R4 $ 
50 FB69 CF44 64FD 0491 91 MULH TABHI Raj, RO 3 Fise MUL by 2**(J/16) 
50 FC62 CF44 60FD 0498 8 ADDH TABLOCLR4], RO 3 Add in LO of 2**(J/16) 
50 FBSB CF44 60FD Rene 7 ADDH TABHILR4], RO 3; And then HI of 2**(J/16) 
56 =osOOF CA 04A6 95 BICL #15, R6 : Is, 
oc 13 Q4A9 96 BEQL ‘ fe I= then don 
56 56 FC 8F 78 OQ4AB 97 ASHL #-4, R6, Rb 5 Position 4 A. addition 
50 56 CO 048 98 ADDL2 = R6, RO ; tg exponent. 
048 99 3 MUL by **] by oupenent addition 
50 B65 0483 400 TSTW RO ; Test for over/underflo 
42 15 0485 401 BLEQ EXCEPT ; See what exception len +f Aang A or = 0 
05 0487 186 20S: RSB ; Otherwise return result in 
0488 40 
0488 404 SMTST: 
14 19 0488 405 BLSS 20$ : Exception if exp+4 > 18 
3F8F BF 54 Bl O4BA 406 CMPW R4, #*X3F8F : Eliminate under flow from APPROX1 
08 19 O4BF 407 BLSS 108 > Bypass if E**ARG = 1 
04C1 408 
04C1 409 ;+ 
04C1 410 ; Use Chebyshev series for small arg 
04C1 411 :- 
04C1 oi¢ 
FO38 CF OF 50 75FD QO4C1 41 POLYH RO,#HXPLN1-1,HXPTB1 :Use Chebyshev series 
04C 414 : last term is 1; this will 
04¢ 415 : give desired overhang. 
05 ot 216 RSB ; Answer is OK, return 
50 08 70FD Q4C9 418 10S: MOH $*#1.0, RO : Ex*X is 1.0, store it 
05 Rete 419 RSB 3; And return 
4CE 420 
4CE 421 
4CE 4 § ; 
Roce ? ? ; Handlers for software detected over/underflow conditions follow 
50 73FD O4CE 425 208:  TSTH RO ; If big ARG > 0 goto OVERFLOW 
2A 18 erg ; os BGEQ OVER 
‘D3 2 g 3 ; Underflow; if user has FU set, signal error. Always return 0.0 
ret 430 UNDER: 
54 DC Q4D 431 MOVPSL R4 ; R4 = user’ ish ron routine’s PSL 
00000000 ' GF 28 43 405 4 ¢ CALLS #0, G*MTHSSJACKET_TST ; RO = TRUE if ons roe. jacket routine 
4 50 § 4DC 4 BLBC : Branch if user did 
54 04 AD 3C O4DF 434 MOVZWL SF$W_SAVE_PSW(FP), R4 =: Get user PSL saved me CALL 


~~ 
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1-002 ATHSHEXP. RO-= Special HEXP peut ine aot 91:32:88 EMTHRTL. SREIATMHEXP MAR: 1 . (5) 
50 7C OQ4E 435 10$: CLRQ RO ; RO = result. LIBSSIGNAL will save in 
4E5 4 : ; CHFSL_MCH_RO/R1 so any handler can fixup 
52 7C OQ4E 4 CLRQ Re D nes 
0D 54 06 €1 OQ4E 4 8 BC #6, R4, 208 ; Has user enabled reese under f Low? 
6—E DD OQO4EB 4 PUSHL ; Yes, return PC from special routine 
7E OO°8F GSA QO4ED 440 MOVZBL #MTHSK_FLOUNDMAT, =-(SP) ; Trap code for hardware floating underflow 
4F1 = 441 ; Convert to MTHS_FLOUNDMAT (32-bit VAX-11 
4F1 rr 3; Exception code) 
00000000 ' GF 02 FB Q4F1 44 CALLS #2, G*MTHS$SIGNAL 3; Signal (condition, PC) 
05 bees ret: 20$: RSB ; Return 
Bere a8 EXCEPT: ; 
56 DS QO4F9 44 TSTL R6 ; Test sign of I; if 1 <0 
D6 =—«19 are re BLSS UNDER : Go to underflow handler 
gene $39 3; Signal floating overflow, return reserved operand, -0.0 
O4FD rd OVER: ; Else process for overflow 
6—E OD O4FD 45 PUSHL (SP) ; Return PC from special routine 
7E OO'B8F QA OQ4FF 454 MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; Hardware floating overflow 
1 OF 79 0503 455 ASHQ #15, #T, RO ; RO/R3 = result = reserved operand 
0507 456 ; 70.0. RO/O will be copied to 
0507 457 ; Signal mechanism vector (CHFSL_MCH_RO/R1) 
0507 458 ; So can be fixed up by any error 
0507 459 3; Handler 
52 7C 0507 460 CLRQ R2 aes 
00000000'GF 02 FB 0509 461 CALLS #2, G*MTHSS$SIGNAL ; Signal (condition, PC) 
05 0510 Pet RSB ; Return - RO restored from CHF$L_MCH_RO/R1 
$811 464 
0511 465 ; End of module MTHSHEXP 
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MTHSHE XP 16-SEP-1984 0 AX/VMS v04-00 Page 
Symbol table at ets rane 3 88 ERTHRTL. Sh SATHH HHEXP.MAR; 1 ° 
ACMASK = 0000407 
APPROX NN R 
EVAL 000477 R 
EXCEPT QOOO4F9 R 
EXP = tt 
HXPLN = 44 F 
HXPLNI = 000 B05 
HXPTAB si 44 00 R 0 
HXPTB1 000 gr R 0 
H_16L0G2_E QOQOOOSFO R 0 
H_LN2_OV_16_HI 00000400 R 0 
H_LN2_OV_16_L0 sated R 0 
ONG = 00000004 
MTHSSJACKET_HND teeeenee 3 
MTHSSJACKET_TST eeeneeee 0 
MTHSSSIGNAL eeeereee 00 
MTHSHE XP 00000420 RG 02 
MTHSHEXP_R6 00000436 RG 02 
MTHSK_FLOOVEMAT eneeere x 00 
shi _FLOUNDMAT tetteeee 00 
OV QOOO04FD R 02 
Oy eu SAVE_PSW = 00000004 
SMTST 00000488 R 02 
TABHI 00000000 R 02 
TABLO 00000100 R Os 
UNDER 00000403 R 0 
x = 00000008 
X_16LOG2E = 00007D12 
err eer ese o se eeeces + 
! Psect synopsis : 
PSECT name Allocation PSECT No. Attributes 
ABS . 00000000 0.) 00 ¢ O.) NOPIC USR CON ABS LCL ao NOEXE NORD —— NOVEC BYTE 
$ 00000000 ( 0.) QO1¢ +1.) NOPIC USR CON ABS LCL NOSHR EXE RD T NOVEC BYTE 
MTHSCODE 00000511 (¢ 1297.) 02 ¢ 2.) PIC USR CON REL LCL SHR’ EXE~ RD NOURT NOVEC LONG 
eeeeeeeeceeeeoeeeeeeeoeoen + 
! : Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 29 00:00:00.1 0:00:00.44 
ponnane processing i :00: §-} + 899-83 
Syabol table sort 0 :00: 5: 6 0:00: 3 
s 2 100 700: ® 228: 4. 
aE table output 4 0 : 8: 8.19 
poset synopsis tee 3 90 88 8 “Be 
ross-reference outpu 00:00. 
Assembler run petals 386 3:0 2:35 8804950 


The working set Limit was 1050 


9258 bytes (19 pages) of virtua 


memory were used to buffer the intermediate code. 
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| VAX-11 Macro Run Statistics 6-SEP-1984 11:24: MTHRTL.SRCIMTHHEXP.MAR; 1 (5 


There were 10 pages of symboi table space rece d to hold 56 non-local and 5 local symbols. 
5 source Lines were read in Pass 1, produci ob object records in Pass 2. 
9 pages of virtual memory were used to define 8 macros. 


) 


trem awne ana na aren moe maws coe n $ 
! ; Macro Library statistics H 
Macro library name Macros defined 
as GETS were required to define 4 macros. 
| There were no errors, warnings or information messages. 


| MACRO/ENABLE =SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=L1S$:MTHHEXP/OBJ=OBJ$:MTHHEXP MSRC$:MTHJACKET/UPDATE=(ENHS$:MTHJACKET) +MSRCS$: 
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